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Abstract 

The microcanonical transfer matrix is used to study the zeros of the partition function of the Q-state Potts 
model. Results are presented for the Yang-Lee zeros of the 3-state model, the Fisher zeros of the 3-state model 
in an external field H q < 0, and the spontaneous magnetization of the 2-state model. In addition, we are able 
to calculate the ground-state entropy of the 3-state model and find sq = 0.43153(3) in excellent agreement with 
the exact value, 0.43152... 



I. Introduction 

The Q-state Potts model p] in two dimensions ex- 
hibits a rich variety of critical behavior and is very fer- 
tile ground for the analytical and numerical investiga- 
tion of first- and second-order phase transitions. With 
the exception of the Q = 2 Potts (Ising) model in the 
absence of an external magnetic field H, exact solutions 
for arbitrary Q are not known. However, some exact re- 
sults have been established for the ferromagnetic Potts 
model Q. For Q < 4 there is a second-order phase tran- 
sition, while for Q > 4 the transition is first order. 

The bond-energy for the Q-state Potts model is (in 
dimensionless units) 



E= £ [l-Sfaaj)], 

<i,j> 



(1) 



where < i,j > indicates a sum over nearest-neighbor 
pairs, crj = 0,...,Q — 1, and E is a positive integer < 
E < Nb, where Nf, is the number of bonds on the lattice. 
In addition to the energy there are Q order parameters 



M„ 



k 



(2) 



where q is a fixed integer between and Q— 1. Note that 
< M q < N s is also an integer and N s is the number of 
sites on the lattice. 

If we denote the density of states with energy E by 
£l(E), then the partition function for the Potts model is 



z(y) 



(3) 



where y = e~@. From (3) it is clear that Z(y) is sim- 
ply a polynomial in y. If we wish to study the parti- 
tion function in an external field which couples to the 
order paramter, (2), then one needs to enumerate the 
states with fixed energy E and fixed order parameter 
M, Q(M,E). The partition function in a magnetic field 
H q is again a polynomial given by 



N s N b 

Z(x,y)= £ £n Q (M,E)x 

M=0 E=Q 



M y E , 



(4) 



where x = e@ Hq . 

II. Microcanonical transfer matrix 

By the microcanonical transfer matrix (//TM)|| it 
is possible to obtain exact integer values for £l(E) and 
Q(M,E). Here we describe briefly the fiTM. on an L x 
L lattice with cylindrical boundary conditions. First, 
an array, uj, which is indexed by the energy E and L 
variables <7j, < i < L — 1 for the first row of sites is 
initialized as 



L-l 



UJ 



(E,a Q ,Gi,...,a L -x) = S(E- K a h <ri+i)), (5) 



i=0 



where h(q,p) is the contribution to the energy from a 
single bond. 

Following the Binder algorithm each spin in the 
row is traced over in turn, introducing a new spin vari- 
able from the next row, 

u:(E;a' Q , ...,ox_i) = ^uj(E - h(a ,a' );ao, ...,0X_i). 

CO 

(6) 

This procedure is repeated until all the spins in the first 
row have been traced over, leaving a new function of 
the L spins in the second row. The horizontal bonds 
connecting the spins in the second row are then taken 
into account by shifting the energy, 

L-l 

uj{E;a' ,... } a' L _ l ) =u(E- ^ Koi,Oi+i);o' ,...,o , L _ 1 ). 

i=0 

(7) 

This procedure is then applied to each row in turn until 
the final row is reached. The density of states is then 
given by 

°(£) = EE- E ^^^..Wl-i)' (8) 



III. Partition function zeros 

The analytic structure of the partition function is 
completely determined by its Yang-Lee (YL) zeros in the 
complex x plane jBI] and its Fisher zeros in the complex y 
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Figure 1: Zeros of a 6 x 6 three-state Potts model in 
the complex x plane for several values of y (cylindrical 
boundary conditions). 



Figure 2: Critical temperatures of the three-state Potts 
ferromagnet function of the magnetic field. The 
upper dotted line is the Ising transition temperature in 



the limit H„ 



-oo, while the lower dotted line shows 



plane f§. Using the ^TM we have studied the YL zeros 
of the two-dimensional Potts model for Q = 2, for which 
the circle theorem of Lee and YangQ applies, and, for 
the first time, for Q > 2^. Figure 1 shows the YL zeros 
for the three-state Potts model. At y = 0.5y c the zeros 
are uniformly distributed close to the unit circle. As 
the temperature is increased the edge singularity moves 
away from the real axis and the zeros detach from the 
unit circle. Finally, as y approaches unity, the zeros site is 111] 
converge on the point x = 1 — Q@. The Q-state model 
in an external field H q < is in the same symmetry 
class as the Q — 1 state model. We have studied the 



the critical temperature of the three-state Potts model 



for H q = 0. 

One of the most interesting properties of this model 
is that for Q > 2 the ground-state is highly degenerate 
and the ground-state entropy is nonzero. However, the 
exact ground-state entropy of the APM is not known 
except for Q = 3, where the ground-state entropy per 



so 



2 3 



0.43152. 



(10) 



Fisher zeros of the 3-state model as a function of 

and determined the critical point and exponents. The 

critical line is shown in Figure 2. 

By introducing the density of zeros g(6,y) Lee and 
Yang[||] showed that the spontaneous magnetization of 
the Ising model is given by 



From the density of states we have obtained the 
exact ground-state entropy for finite lattices. Table 1 
shows the ground-state degeneracy Qq(L) of the L x 
L three-state APM. The ground-state entropy is then 
given by 

hiQ (L) 



so(L) 



(11) 



m (y) = 2vrc/(0,y). 



(9) 



Even though they have calculated the density of zeros 
of the one-dimensional Ising model and confirmed the 
well-known fact that the model does not exhibit ferro- 
magnetism, the exact form of the density of zeros is not 
known for d > 2. We have studied the density of zeros 
of the two-dimensional Ising model numerically^, and 
using (9) we have calculated the spontaneous magneti- 
zation (Figure 3) which agrees very well with the exact 
result g. 

IV. Ground-state entropy 

The antiferromagnetic Potts model (APM) is much 
less well understood than its ferromagnetic counterpart. 



By using the Bulirsch-Stoer (BST) algorithm [|i~2jl , we 
have extrapolated our results for finite lattices to infinite 
size. The BST estimate for Q = 3 is so = 0.43153(3), in 
excellent agreement with (10). 
V. Discussion 

Using the ^TM, we have calculated the restricted 
density of states £l(M, E) from which the partition func- 
tion at any temperature and in an arbitrary magnetic 
field can be evaluated. We have studied the Yang-Lee 
zeros of the 3-state model and demonstrated that, un- 
like the Ising model, the locus of zeros is not a circle. 
Following Lee and Yang we have used the density of ze- 
ros to calculate the spontaneous magnetization for the 
2-state model and we are currently extending this work 
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Figure 3: Extrapolated density of zeros as a function of 
y. The solid curve is Yang's exact result for the spon- 
taneous magnetization. Open circles are our results for 
cylindrical boundary conditions, and open squares are 
for periodic boundary conditions. 

Table 1: The ground-state degeneracy of the three-state 
APM for even-size lattices (cylindrical boundary condi- 
tions). 



L 



n (L) 



4 4626 

6 37284186 

8 9527634436194 

10 77048019386429374638 

12 19698820973096872077077373450 

14 159147870862104841838351532192943853490 

to Q > 2. Since the zeros of the partition function en- 
compass both the ferromagnetic and antiferromagnetic 
regimes, we are currently investigating the critical be- 
havior of the antiferromagnetic Potts model. Also, the 
/iTM is easily extended to three-dimensional lattices, 
and we are now extending all of the above calculations 
to cubic lattices. 

Recently it has been suggested |b| that the micro- 
canonical entropy may exhibit weaker finite-size effects 
than the free energy. The microcanonical ensemble also 
facilitates the study of phase coexistence and the study 
of phase transformations in cluster and nuclear physics fll3| 
The //TM method yields exact values for the micro- 
canonical entropy of finite systems allowing us to study 
finite size effects in the microcanonical ensemble. 



[2 



[10 



[11 



R. B. Potts, Proc. Camb. Phil. Soc. 48 (1952) 106; 
F. Y. Wu, Rev. Mod. Phys. 54 (1982) 235; R. 
J. Creswick and S.-Y. Kim, |cond-mat/970101~8 ; J. 
Phys. A 30 (1997) 8785, and references therein. 

L. Onsager, Phys. Rev. 65 (1944) 117; C. N. Yang, 
ibid. 85 (1952) 808. 

B. Stosic, S. Milosevic, and H. E. Stanley, Phys. 
Rev. B 41 (1990) 11466; G. Bhanot, M. Creutz, 
and J. Lacki, Phys. Rev. Lett. 69 (1992) 1841; R. 
J. Creswick, Phys. Rev. E 52 (1995) 5735. 

K. Binder, Physica 62 (1972) 508. 

C. N. Yang and T. D. Lee, Phys. Rev. 87 (1952) 
404; T. D. Lee and C. N. Yang, ibid. 410. 

M. E. Fisher, in: Lectures in Theoretical Physics, 
vol. 7c, W. E. Brittin, ed. (University of Colorado 
Press, Boulder, 1965) p. 1. 

S.-Y. Kim and R. J. Creswick, Phys. Rev. Lett. 81 
(1998) 2000. 



S.-Y. Kim and R. J. Creswick, |cond-mat/9806294 
(Phys. Rev. E, in press). 

R. J. Creswick and S.-Y. Kim, Phys. Rev. E 56 
(1997) 2418. 

R. Shrock and S.-H. Tsai, Phys. Rev. E 55 (1997) 
5165; S.-H. Tsai, ibid. 57 (1998) 2686, references 
therein. 

E. Lieb, Phys. Rev. 162 (1967) 162. 



[12] R. Bulirsch and J. Stoer, Numer. Math. 6 (1964) 
413; M. Henkel and G. Schiitz, J. Phys. A 21 (1988) 
2617. 

[13] A. Hurler, Z. Phys. B 95 (1994) 63; M. Kastner, M. 
Promberger, and A. Hurler, in: Computer Simula- 
tion Studies in Condensed-Matter Physics, vol. 11, 
D. P. Landau and H.-B. Schiittler, eds. (Springer, 
Berlin, in press); D. H. E. Gross, Phys. Rep. 279 
(1997) 119, and references therein. 



3 



